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I/-) ■ ABSTRACT 

Aims. New data and results on the optical behaviour of the prominent blazar PKS 0735+178 (also know as OI 158, S3 

■ 0735+17, DA 237, 1ES 0735+178, 3EG J0737+1721) are presented, through the most continuous BVRI data available 
^ ■ in the period 1994-2004 (about 500 nights of observations). In addition, the whole historical light curve, and a new 

' photometric calibration of comparison stars in the field of this source are reported. 

£\J , Methods. Several methods for timeseries analysis of sparse data sets are developed, adapted and applied to the recon- 

, structed historical light curve and to each observing season of our unpublished optical database on PKS 0735+178. 

• Optical spectral indexes are calculated from the multi-band observations and studied on long-term (years) durations 

as well. For the first time in this source, variability modes, characteristic timescales and the signal power spectrum are 
£ — ( explored and identified over 3 decades in time with a sufficient statistics. The novel investigation of mid-term optical 

scales (days, weeks), could be also applied and compared to blazar gamma-ray light curves that will be provided, on 
the same timescales, by the forthcoming GLAST observatory. 
i-Q ' Results. In the last 10 years the optical emission of PKS 0735+178 exhibited a rather achromatic behaviour and a vari- 

ability mode resembling the shot-noise. The source was in a intermediate or low brightness level, showing a mild flaring 
activity and a superimposition/succession of rapid and slower flares, without extraordinary and isolated outbursts but, 
at any rate, characterized by one major active phase in 2001. Several mid-term scales of variability were found, the 
more common falling into duration intervals of about 27-28 days, 50-56 days and 76-79 days. Rapid variability in the 
historical light curve appears to be modulated by a general, slower and rather oscillating temporal trend, where typical 
amplitudes of about 4.5, 8.5 and 11-13 years can be identified. This spectral and temporal analysis, accompanying our 
data publication, suggests the occurrence of distinctive signatures at mid-term durations that can likely be of transitory 
JJ^j ' nature. On the other hand the possible pseudo-cyclical or multi-component modulations at long times could be more 

stable, recurrent and correlated to the bimodal radio flux behaviour and the twisted radio structure observed by several 
years in this blazar. 
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1. Introduction 

The rapid and violent optical variability is one of the 
defining properties of blazars, and variability studies are 
important in understanding the physics of AGN in gen- 
eral. Characteristic timescales, fluctuations modes, flares 
shapes and amplitudes, duty cycles and spectral changes, 
correlations and temporal lags between variations in dif- 
ferent spectral bands, provide crucial information on the 
nature, structure and location of the emission compo- 
nents and on their interdependencies. In particular the so 
called low/intermediate-frequency peaked BL Lac objects 
(LBL/IBL), have the peak of the synchrotron emission 
around infrared and optical wavelengths and commonly 



Send offprint requests to: S. Ciprini 
e-mail: stefano.ciprini@utu.fi 



show large-amplitude flares characterized by prominent flux 
variations in a wide range of temporal scales. The rapid 
optical variations of LBL and IBL are also systematically 
larger and with shorter duty cycles than those of the high 
energy peaked BL Lac objects (HBL). Hence a multi-band, 
possibly well sampled and extended optical monitoring is 
an important and subsidiary element of the standard multi- 
wavelength (MW) analysis. MW observing campaigns pro- 
vide, more or less, short snapshots of the targets, lacking of 
information about their mid/long-term evolution. Even if 
the optical band has a narrow spectral extension, it can 
yield useful information about the synchrotron emission 
peak and possible disk/host-galaxy contributions. Moreover 
long-term (historical) records of blazar variability are avail- 
able at optical wavelengths for several bright objects, al- 
though data collected in the past are rather sparse. Small- 
size and dedicated (possibly automatic) telescopes, in con- 
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junction with international consortiums, have recently in- 
creased the amount of photometric data, sometimes with 
a fair continuous sampling during specific observing cam- 
paigns. 

In this paper we present more than 10 years (Feb. 1993- 
Feb.2004) monitoring data about the blazar PKS 0735+178 
(1332 photometric points in four BVRI Johnson-Cousins 
filters, obtained during about 500 observing nights). Our 
effort represents the best optical monitoring available for 
this object regarding to continuous and long-term coverage. 
This optical programme allowed to study colours and the 
continuum spectrum (mainly in VRI bands) and, for the 
first time, enabled to study mid-term scales (days, weeks), 
over an extended data set. These timescales were almost un- 
explored in blazars due to the irregular/poor sampling and 
the low statistics in the optical regime. Time series anal- 
ysis accompanying our data publication, is performed for 
both our observations and the historical light curve (1906- 
2004), while a new photometric calibration of comparison 
stars in the field of the source is also reported, as useful 
reference for future optical observations and monitoring. In 
the historical light curve there are obvious differences in 
data quality, accuracy and sampling over time, that can 
yield biases, noise, spurious and fakes signatures. However 
we remark that the last 33 years portion (1970-2004) of 
the historical light curve holds a sufficiently regular sam- 
pling to allow meaningful statistical results on long-term 
intervals too. Data binning when needed, the employment 
and comparison of 7 different temporal analysis methods 
suitable for unevenly sampled data set, and the calcula- 
tion of the power spectrum given by the gaps, secure us 
to have determined and reported only real and intrinsic 
time signatures. We note finally that the main aim of our 
paper was to investigate the variability behaviour on such 
intermediate scales through our 10-year observations. Data 
published in this work were obtained by 4 optical obser- 
vatories: the Perugia University Observatory (Italy), the 
INAF Torino Observatory (Italy), the Tuorla Observatory 
(Finland), and the Sabadell Observatory (Spain). Perugia, 
Torino and Sabadell data on PKS 0735+178 are unpub- 
lis hed, while part of t he Tu orla data were alrea dy published 



BL Lac object in ICarswell et al.l ([1974I ). This s ource is op- 
ticall y bright, highly vari able, and both ra dio (jKiihr et al.l 
I1981D and X-ray selected ([Elvis et al.lll992ft . PKS 0735+178 
has been extensively studied in the radio regime. The ra- 
dio flux appear to vary quite slowly wit h some outbursts 



([Baath fe Zhang! 119911: 

1999: iTerasranta et al 



Terasranta et al.lll992t lAller et alj 
12004), 



but the re is not any 
evide n ce of a radio-optical corre l ation (|Clements et al.l 
119951: lHanski. Takalo. fe Valtaoial 12002!) or periodicity 
(|Ciaramella et al.ll2004h . Early radio observations of PKS 
0735+178 showed a peculiar spectrum soon interpreted 
as the superposition of incoherent synchrotron radiation 
emitted by distinct and homogeneous radio components, 
conspiring to add up to an overall very flat shape (the 
source was indeed nicknamed as the "Co smic Conspiracy" 
lMarscherlfl977l 119801 : ICotton et al.lll980h . Several moving 
components and an unusual, complex morphology charac- 
terized by a twiste d jet were observed in VLBA/VLBI radio 
imag i ng (see, e.g. iPerlman fe Stockd 11994 iGabuzda et alJ 



1994lGabuzda. Gomez, fe Agudd|200lt lG6mez et al.ll200lb 



20021: lOihaet al 



2004 



iKellermann et all 



Kataiainen et al.l ((2000). Optical data from lQian fe Tad 



Ho man et al. 

I2OO4J). PKS 0735+178 has one of the most bent radio 
jets among AGN observed at milliarcsecond (mas) scales. 
A bimodal scenario in which periods of enhanced activ- 
ity with ejection of superluminal components are followed 
by epochs of low activi ty with a highly twisted jet ge- 
omet ry was suggested ([Gomez et al.1 l200lt lAgudo et al.l 
l2006h . PKS 0735+178 is also a X-ray and gamma-ray 
(EGRET, 3EG J0737+172 1) emitting blazar dKubo et al.l 
119981 : lHartman et alj [i~999l see, e.g. ). The spectral en- 
ergy distribution (SED) evinced PKS 0735+178 as a low 
or intermediate-energy peaked BL Lac object (LBL/IBL), 
where the IR-UV synchrotron continuum dominates the to- 
tal observed power. The very low X-ray variability with 
respect to the high optical-IR variations, supported the 
idea that X-rays are produced by inv erse Compton mecha- 
nism in some mas radio com ponents (jBregman et al.lll984 
iMadeiski fe Schwartz! [l988). The 7-ray flux of this blaza r 
appeared likewise no strongly variable (jNolan et al.l l2003). 

The optical spectrum of PKS 0735+178 shows an ab- 
sorption line due to an intervening system at 3980 A that, 
if identified wit h Mg-II, provides a lower redshift limit 



(|2004) have been also added to improve a few the sampling. f z > 0,424 (ICarswell et all 11974 iBurbidge fe Hewitt! 



The paper is organized as follows: in Sect. [2] we review 
briefly the optical knowledge about PKS 0735+178, while 
in Sect. [3] we mention the observing and data reduction 
techniques. A new photometric calibration of comparison 
stars in the field of the source is presented in Sect. HI and 
the BVRI light curves collected during our monitoring are 
showed in Sect. [5] The reconstructed historical light curve 
is described in Sect. [H] while in Sect. [7] the analysis of the 
multi-band behaviour is reported computing the optical 
spectral indexes. A joint temporal analysis of our data and 
the historical light curve is performed in Sect. [HI and sum- 
mary and conclusions are outlined in Sect.[9l 

2. Optical properties of PKS 0735+178 

The radio object PKS 0735+178, belonging to the Parkes 
catalog (other most used names are: PKS 0735+17, S3 
0735+17, OI 158, DA 237, VRO 17.07.02, PG 0735+17, 
RGB J0738+177, Uy 0735+17, RX J0738. 1+1742, 3EG 
J0 737+1721) w as identified with an optical point source 
bv iBlakd (|l970h . Afterwards it was classified as a classical 



19871 : iFalomo fc Ulrichl 120001 : iRector fc Stockd l2001h . A 



strong Lyman- alpha absorption line has also been de- 
tected by the IUE satellite at t he same redshift 
(jBregman. Glassgold. fc Hugginsl Il98lh . This absorption 
was not identified in deep optical imaging, even if a 
very faint emission was detected about 3.0-3.5" NE/E 
(projected distance 22-25 kpc a t z = 0.424 ) from 
the object (jFalomo fc Ulrichl l2000fc iPursimo et all I2002D . 
The host galaxy of PKS 0735+178 remains unresolved 
in optical imaging ([Scarpa et al.1 120001 : IFalomo fc Ulrichl 
I2OOO1 : IPursimo et al.l 20021 ). but the source has two well 
resolved companion galaxies. The galaxy at 7" NW, 
declared distorted by interaction with PKS 0735+178 
(|Hutchings. Johnson, fc Pvkdll988l ). does not show marks 
of interaction in more recent and higher resolution 
images, and a redshift of z = .645 was obtained 
for it (IStickel. Fried, fc Kuehrl 119931: IScarpa et all 120001 : 
IFalomo fc Ulrichl l2000h . In addition a brighter galaxy lo- 
cated at 8.1" SE cannot be the absorber due to the 
great projected distance from our blazar. A third, very 
faint, elongated structure at 3-3.5" NE was detected as 
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periodica l components were claimed, with values of 1.2 , 
4.e_ycars (ISmith. Leacock. fc Webb l l 19881: 1 Webb et al.fl 988: 



Fig. 1. Finding chart with the optical comparison stars 
identified in a field of 15'xl5' centered on PKS 0735+178 
(within the double bar). The field is elaborated from a 
frame of the Digitized Sky Survey. The new photometric 
VRI calibration of the stars CI, C, D, C2, C4, C7, and A 
is reported in Tablc[T] A, C, and D stars belong a lso to the 
photometric sequence calibrated bv lSmith et a .1 (Il985h . to 
the photometric sequence published in IWind (|l973l . here 
with a B sequ e nce t oo), and the sequence appeared in 
IVeron fe Veronl (|1975L A and D stars only). Stars C2, C4, 
C7 belong also t o the p hotometric sequence calibrated by 
iMcGimsev et al.l (|l976l . named there 1,2,3), while for the 
CI star the calibration is totally new. 



well (|Pursimo et al.lll999tlFalomo fc Ulrichll2000D . and this 
could be related to the intervening absorption at z = 0.424. 
The lower limit for the redshift z > 0.5 obtained assu ming 
typical properties for the host (|Falomo fc UlricrJl2000h is 
consistent with the limit derived from intervening absorp- 
tion. 

When combined with our data the historical optical 
curve of PKS 0735+178, s tarting from JD 2417233, i.e. Jan. 
22, 1906 (|Fan et alJll997h . spans almost 100 years. Optical 
vari ations are often o f larger amplitudes than the infrared 
one (jFan fc Lirj|2000t ) . Correlations between t he spectral in- 
dex and the optical brigh tness were observed (jSitko fc Sitkd 
Il99lt iLin fc Fan] Il998f ). alongside a spectral flattening 
(blueing) with t he source brightening (jBrown et al.lll989t 
ILin fc Fan1ll998h . On the other han d this type of correlation 
appear to be weak (iGu et al.ll2006h or oppo site (i.e. spectral 
steepening, reddening; iGhosh et all [2000) in other multi- 
band observations. The large st optical variations registered 
are o n the order of 3-4 mag ([Pollock et alJll979l : I Fan et aP 
1 1997T) at long (years) ranges. Some intra-day (IDV) and 
inter-day variations until 0.5 mag were repor t ed in the 
optical history of this blazar dXie et al.l fl992i : I Fan et aP 
Il997t iMassaro et aPll995l : IZhang et al.ll2004[ ). Over the pe- 
riod 1995-1997, the optical IDV a nd microvariatio ns were 
rare and with a small amplitude (|Bai et al.l Il999h . while 
no clear evidence was found in more recent observations 
(jSagar et al.l [2004). Several possible recurrent and pure 



ISmith fc Nairlll995l) : 14.2, 28.7 yearsfean et al.lll997h : 8.6, 
13.8, 19.8, 37.8 years (jQian fc Tadl2004h . evenif we remark 
that such scales are derived by different data sets and dif- 
ferent epochs. The long-term analysis performe d with the 
Jurkeyich's method on a more complete data set (|Fan et al.l 
Il997t iQian fc Tao 2004) postulated a main periodical com- 
ponent of about 13.8-14.2 years, but our temporal analysis 
(Section^]) suggest other and shorter long-term signatures. 

PKS 0735+178 has also a relatively high degree of op- 
tical polarization showing very different levels covering the 
whol e range f rom a bout 1% up to 30% (see, e.g. Mead et all 
199 rtlTakaldll991l:lTakalo et al]ll992t IValtaoia et al.lll99lL 



1993 ; Tommasi et al.ll200l[) . Only a modest variability of 



this optical polarization was observed on inter/intra night 
durations, that was interpreted as owed to substructures 
of different polarization and variable intens ity in the jet. A 
prefe rred polarization level over few years (jTommasi et alj 
l200lh could indicate quiescence and stability in the under- 
lying jet structure. 



3. Observations and data reduction 

Photometric observations were carried out with four tele- 
scopes. The Newtonian f/5, 0.4 m, Automatic Imaging 
Telescope (AIT) of the Perugia University ObservatorjQ, 
Italy (451 meters above sea level, a.s.L), a robotic 
telescope equipped with a 192 x 165 pixels CCD 
array, thermoelectrically cooled with Peltier elements 
(|Tosti. Pascolini. fc Fiorucclll99fih . The REOSC f/10, 1.05 
m, astrometric reflector of the Torino Observatorj0, Italy 
(622 meters a.s.L), mounting a 1242 x 1152 pixel CCD ar- 
ray, cooled with liquid nitrogen and giving an image scale of 
0.467" per pixel. The Dall-Kirkham f/8.45, 1.03 m reflector 
of the Turku University Tuorla ObservatorjQ, Finland (60 
meter a.s.L), equipped with a 1530 x 1020 pixel CCD cam- 
era, thermoelectrically cooled. The Newtonian 0.5m tele- 
scope of the Sabadell ObservatorjQ, Spain used in two inter- 
changeable configurations (Newton at f/4, and Cassegrain- 
Relay at f/15), and equipped with a 512 x 512 FLI CM-9 
CCD array. The Perugia and Torino telescopes were pro- 
vided with standar d BV (Johnson) and R r Jr (Cousins) fi l- 
ters (<Bessel]| 119791 iFiorucci fc Munaril [20031 : iBesseill 12001 . 
The Tuorla lm telescope, and the Sabadell telescope were 
equipped with V and R c filters. 

All the observatories took CCD frames and performed 
a first automatic data reduction using standard methods, 
to correct each raw image (for dark and bias background 
signals where needed) and to flat fielding, to recognize the 
field stars, and to derive instrumental magnitudes via aper- 
ture photometry (Perugia and Tuorla) or circular Gaussian 
fitting (Torino). The single frames are then inspected to 
evaluate the quality of the image, the reliability of the 
data, and to search for spurious interferences. Comparison 
among data obtained with these different telescopes on the 
same night reveals a good agreement, and no detectable 
offset is found. The matching with dat a taken in simulta- 
neous epochs at Shanghai Observatory (Qi an fc Tad [2004) 
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Table 1. A new optical photometric calibrations for field of PKS 0735+178, with VR C I C Johnson-Cousins sequences of 
comparison stars CI, C, D, C2, C4, C7 and A (see Fig. Q]). Data are obtained at the Perug ia University Obse r vatory 
and a re adopted in this work . The previous avai l able c o mparison sequence s were published in ISmith et al.l (|l985f) ; [Wing 
(|!973h : IVeron fc Veronl (fl975h : iMcGimsev et all (fl976l ): ISmith et all (|l985h . Magnitudes of common stars are roughly in 
agreement within the uncertainties. The VR C I C sequence reported in this table, joint with the U, B values reported by 
ISmith et alJ (11985D . is suggested for future optical observations and monitoring. 



PHOTOMETRIC SEQUENCES FOR. PKS 0735+178 COMPARISON STARS 



Star R.A. Dec. V R c I c 

(J2000.0) (J2000.0) [mag] [mag] [mag] 

CI 07 38 00.5 +17 41 19.9 13.26 ± 0.04 12.89 ± 0.04 12.57 ± 0.04 

C 07 38 02.4 +17 41 22.2 14.45 ± 0.04 13.85 ± 0.04 13.32 ± 0.04 

D 07 38 08.3 +17 44 59.7 15.90 ± 0.05 15.49 ± 0.05 15.12 ± 0.06 

C2 07 38 08.5 +17 40 29.2 13.31 ± 0.04 12.79 ± 0.04 12.32 ± 0.04 

C4 07 38 11.6 +17 40 04.4 14.17 ± 0.05 13.80 ± 0.04 13.48 ± 0.04 

C7 07 38 20.7 +17 40 51.2 15.01 ± 0.06 14.70 ± 0.06 14.37 ± 0.05 

A 07 38 23.4 +17 42 43.0 13.40 ± 0.05 13.10 ± 0.05 12.82 ± 0.05 



showed a good agreement too (see Fig. [3]). The precision 
level in the light curve of PKS 0735 assembled in this way 
is enough for a variability analysis performed on interme- 
diate and long-term timescales. Moreover the time series 
analysis of shorter timescales (< 200 days) is performed in 
each single observing season using i?-band data, that were 
mainly obtained by 1 telescope (see Tab.[2|). 

4. Comparison stars photometry 

Calculation of the source magnitude is easily obtained by 
differential photometry with respect to comparison stars in 
the same field of the object. The discussion of the adopted 
comparison star sequence is crucial for analysis of opti- 
cal data obtained during blazar monitoring observations, 
as the photometric sequence affects data quality and re- 
liability. In order to obtain a dependable photometric se- 
quence for PKS 0735+178, we selected a set of non-variable 
stars with brightness comparable to the object and differ- 
ent colours (see the finding chart in Fig. [T]) . Photometric 
calibrations of these stars were derived from 13 optimal 
photometric nights between 1994 and 1996 at the Perugia 
University Observatory using Landolt standards. The sta- 
bility of the sequence for the stars CI, C, D, C2, C4, was 
well tested and verified during data reduction of the over- 
all database (for stars A and C7 there are less measure- 
ments because of the small FOV of the instrument). Our 
new photometric sequence is presented in Table [T] showing 
the V (Johnson), and R C ,I C (Cousins) photometric values. 
Pr evious calibr at ions of this sta r -field were performedin 



I Wind (fl973h: IVeron fc Veronl (|1975h : iMcGimsev et~aTl 
(|!976h : ISmith et all (119851). In " particular the determina- 
tions of IMcGimsev et al.l (|l976l ) are photoelectric in the 
UBV Johnson system, thus well similar to our calibrations 
in B, V bands using a CC D detect o r. Dis crepancies are 
small if the spec ifications o f Bessell (1990) are respected 
and determined (|Fiorucci fc Munarif 20031 ). 

This new photometric calibration of the PKS 0735+178 
field, (suitable also for telescopes with small FOV), is 
slightly more extended and accurate with respect to the 
past calibrations. Table [T] reports the V, R c , I c values for 
seven comparison stars, while the sequence for the star 
denoted with CI is completely new. In each photometric 
night standard Landolt stars were observed at different air- 



masses, and the calibration line as a function of the airmass 
was constructed (neglecting the colour corrections being al- 
ways smaller than the instrumental errors). The standard 
magnitudes of such stars were derived (with error equal to 
the quadrature sum of the linear regression error and the 
instrumental error on the single star). The typical error 
for each night and each star is between 0.03 and 0.1 mag 
(depending on the luminosity of the star and the atmo- 
spheric conditions) . Data shown in Table [T] are the result 
of weighted averages on the values of each night (weight 
equal to 1/c 2 ), whereas the error estimation is equal to the 
standard deviation weighted on the averages. This uncer- 
tain is higher than the standard deviation on each single 
night, because of some systematic errors different on each 
night. We chose to report a reliable calibration with reliable 
errors in Table [T] with respect to more accurate but more 
doubtful smoothed values. Colour transformations to com- 
parison stars were not applied, because from the analysis of 
Landolt stars it was not possible to separate this effect from 
the instrumental statistical errors given mainly by the effec- 
tive limits of the Perugia instrument and site. Then again 
the photometric system of the Perugia telescope was devel- 
oped to follow at best the standard Johnson-Cousins sys- 
tem. Finally we note that some of the stars listed in Table 
[Dare quite red (V — I index ranges from +0.58 to +1.13 
mag, suggesting that they have spectral types F, G, possi- 
bly K), but these colour indexes are similar to the colour 
indexes of PKS 0735+178 (average V — I = +0.95 ± 0.1, 
average B — R = +0.91 ± 0.09) hence colour effects arc 
correspondent, and -B-band data were obtained only with 
the larger lm Torino telescope. A detailed description of 
the observing and data reduction procedures, filter system, 
software adopted in the ca librations, and comparison wit h 
other works can be foun d in lFiorucciT Tosti. fc Riz"^il(|l998D ; 
iFiorucci fc Tostil (|l996l ). 

5. Optical light curves from 1993 to 2004 

We monitored the BL Lac object PKS 0735+178 in the 
four B,V,R C ,I C optical bands for more than 10 years, 
from February 2, 1993 to February 17, 2004 (JD=2449021- 
2453053). A total of 1332 BVRI reduced and validated 
photometric points were obtained over a period of 4032 days 
(see Fig. [2]). In order to obtain a more complete light curve, 
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Fig. 2. BVRI magnitude light curves of PKS 0735+178 from 1993 to begin ning of 2004. Dat a cames from our 10- 
year observing monitoring. Published observations from Shanghai Observatory (jQian fc Taoll2004D are added in order to 
improve the sampling. Data sets of different observatories are in agreement within the uncertainties. 



data from the Shangha i Obse rvatory (Jan. 1995-Dec. 2001) 
are added lOian &: Taol (|2004t ). In the best sampled band 
(the i?-band, analyzed in detail in Sect [8]), 709 photomet- 
ric points were collected over 12 observing seasons (see 
Tab[5] and Tab|3]) with 459 nights in total having at least 
one R data point. The last and best sampled 10 observing 
seasons (from the III to the XII, i.e. from October 1994 
to February 2004) have a duration spanning from 144 to 
203 days Tab[3|), an average number of data points per 
night equal to 1.5, an average empty gap between subse- 
quent observations of 3 days, and an average coverage of 
nights with data respect to each season duration of about 
27%. Practically speaking such numbers mean that data 
are not clustered or bunched, and that an enough regu- 
lar monitoring was performed (when permitted by atmo- 



spheric/technical conditions). The priority of our observing 
programme during these years was to perform a constant 
and possibly uniform optical monitoring. Consequently, for 
the first time in PKS 0735+178, this allowed to obtain 
data suitable for a deep and detailed statistical analysis on 
days/weeks timescales. The majority of such i?-filter ob- 
servations were obtained by only two telescopes (over the 
69% obtained by the Perugia telescope, and a further 21% 
by the Torino telescope) having known inter-instrumental 
offsets below 0.1 mag in this band (e.g. IVillata et al.| [2002; 
iBottcher et al.ll2005h . Three examples of such R-band sea- 
sonal light curves with the accompanying time series anal- 
ysis functions are reported in FigfTUl FiglTTI and FigfT2"l 

A direct visual inspection of our 10-year multiband 
light curve (Fig. [2]), show an average optical brightness 
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Fig. 3. A portion of the V-mag light curve where is possible 
to see in more detail data from the different observatories. 

Table 2. The number of photometric BVRI data points of 
PKS 0735+178 obtained by each observatory in the period 
1993-2004 and published in this paper. In the bottom panel 
a summary of the basic information and statistics about our 
data in each band. 

DATA POINTS PER OBSERVATORY 



Obs. 


B 


V 


R 


I 


Tot. 


Period 


Perugia 





226 


490 


281 


997 


Fcbl993-Fcb2004 


Torino 


75 


38 


150 





263 


Decl994-Apr2002 


Tuorla 





55 








55 


Octl995-Feb2001 


Sabadell 








17 





17 


Dec2001-Feb2004 


Shanghai 





115 


52 


138 


305 


Janl995-Dcc2001 


Total 


75 


434 


709 


419 


1637 





1993-2004 DATA STATISTICS 



outburst (R > 14 in the whole data set). Luminosity 
drops/increasing of about 2 magnitudes were common on 
time intervals smaller than half year. From the end of 
1997, a slow increase of the average brightness was clearly 
detected (e.g. the i?-band magnitude never dropped to 
values higher than 16 from beginning of 1998), while in 
2001 a clear brightening phase can be well identified. This 
moderate-outburst phase can be considered comparable to 
the other outbursts seen in the optical history of PKS 
0735+178 (see, Sec. [SJ . 



6. The historical light curve 

The optical history of PKS 0735+178 together with our 
data (Fig. [J) extends over almost one century (from Jan 
22, 1906, JD 2417233, to Feb 17, 2004, JD 2453053). The 
older points in the light curve were obtained using plates 
(mainly from the Landessternwarte Heidelberg-Konigstuhl 
Observatory, G ermany, and the Rosemary Hill Obser vatory, 
Florida, USA IZekl et all fT98lt IWebb et all 11988ft . from 
which a photographic magnitude m pg can be extracted 
and converted in the photometric B magnitude following a 
semi-empirical correction (see, e.g. lLulll972t lKidgerl ll989L 
More recent data have been obtained directly with photo- 
electric or CCD instrum ents. The historica l data collection 
was taken directly from (jQian fe Taoll2004ft with few addi- 
tions and appending our original and derived B-magnitudes 
(the derived B-band data are estimated from our best sam- 
pled i?-mag data after 1993, using a constant colour in- 
dex with value equal to previous works B — R = 0.993, 
iFan et all fl997t iQian fc Tad 120041 ). In general a pruden- 
tial error estimation (taking into account different offsets, 
different data quality, systematic errors and instrumental 
dispersion), needs to be figured out, before to use heteroge- 
neous historical optical light curves for a quantitative anal- 
ysis. In this specific case the further errors introduced by 
rising a constant conversion index from the R and B band 



B 



V 



R 



I 



Total data points 75 434 709 419 

Start date [JD-2449000] 698 45 21 420 

End date [JD-2449000] 3354 4053 4053 4053 

Total period JV to t [days] 2657 4001 4032 3633 

Nights with data N on 52 297 459 259 

N on /N to t fraction 0.019 0.074 0.171 0.071 

Mean num. points x night 1.44 1.46 1.51 1.62 

Total mean gap At [days] 35.9 9.3 5.8 8.7 

Longest gap [days] 780 352 375 356 

Average brightness [mag] 16.319 15.760 15.301 14.693 

Max brightness [mag] 15.863 14.544 14.16 13.59 

Min brightness [mag] 17.453 16.94 16.87 15.97 

Variab. range Am [mag] 1.59 2.39 2.71 2.38 

Absorption coeff. 1 " [mag] 0.152 0.117 0.094 0.068 

Data standard deviation 0.256 0.368 0.515 0.453 

Data skewness 1.23 0.386 0.329 0.155 

Data kurtosis 3.791 1.087 0.019 0.400 

Max flux [mJy] 2.21 6.1 7.3 9.9 

Min flux [mJy] 0.51 0.67 0.60 1.1 



-has also to be counted on. This estimation is difficult with- 
out all the original data sets (plates, frames, etc.), but bas- 
ing on experience a reasonable and prudential upper limit 
to the errors in Fig0] might be considered around the value 
of ±0.4 mag. 

The largest outbursts or brightening phases (mag B < 
15) occurred in the period Dec.1937-Feb.1941, Apr. 1949- 
Feb.1950, around Feb.-Mayl977, in the period Oct. 1980- 
Mar.1981, and Feb.2001-Oct.2001. The brightest outburst 
was observed around mid of May 1977 (JD 2443277-78), 
when PKS 0735+178 reached its historical optical maxi- 
mum (B = 13.9). A rather humped, swinging and oscillat- 
ing long-term trend appear to modulate the rapid variabil- 
ity of PKS 0735+178. This aspect might suggest a cycli- 
cal or intermittent trend, with possibly pseudo-periodic or 
multi-component oscillations running on long-term ranges. 



=7. Optical spectral indexes 



| Values for the galactic extinctio n by NED database 
ijSchleeel. Finkbeiner. fc David 11998ft . 



placed at a mid or low levels, displaying rapid variability 
with a moderate flaring but no extraordinary big/isolated 



The continuum spectral flux distribution of blazars in 
the optical range can be analyzed to distinguish prop- 
erly all the emission components that, together with syn- 
chrotron radiation, contribute to the observed spectrum 
shape. Moreover optical flux variations in blazars are fre- 
quently associated by changes in the spectral shape. This 
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Fig. 4. The historical optical li ght curve of PKS 0735+178 in B band reconstructed by published data (data points 
mainly from the compilation of Qian fc Taoll200"4l . with more few data from other literature) and adding our original 
B and i?-derived data (see text). Error bars are not represented for clarity. The total historical time series obtained in 
this way is composed of 1725 final data points. At least five main outbursts followed by a general humped and rather 
oscillatory long-term trend can be visually recognized. This behaviour is identifiable better in the 1970-2004 portion 
(lower panel), thanks to a higher sampling (lower panel). 



can be revealed by analyzing the magnitude color indexes 
or the flux spectral indexes. In calculating the color in- 
dexes and the continuum spectral slopes, we selected the 
more accurate multi-band data (3 filter at least) about 
PKS 0735+178, obtained by a single telescope, and cou- 
pling frames with a maximum time lag of 20 minutes 
(in order to reduce possibly intrinsic/extrinsic, instrumen- 
tal micro-variations). Since both the host galaxy of PKS 
0735+178 and the feature possibly responsible for the in- 
tervening absorption at z = 0.424 are rather faint (see 
sect. [2]), it is reasonable to neglect the galaxy color in- 
terference and any thermal contribution in the observed 
continuum optical spectra. The observed magnitudes were 
transformed into flux densities, corrected by the Galactic 
absorption (derived by ISchlegel. Finkbeiner. fc Davis!l 998. 
see Tab. [2] lower panel) for the source (located at moder- 
ate Galactic latitude, b=18.07). The absorption is rather 
small (i.e. B — I = 0.8 mag), accordingly the colour correc- 
tion is little in comparison to the mean B — I value of the 
source. Fluxes relative to zero-magnitude values are taken 



from the Johnson- C ousi ns system calibration pres ented in 
iBessefll (|1979L Il990t ) and lFiorucci fc Munaril (|2003f ). 

The optical spectral energy distribution (SED), can be 
expressed conveniently by a power law vF v oc v~ a+ , {y 
being the frequency of radiation and a the spectral in- 
dex). In the optical regime the degree of correlation be- 
tween a and the flux sheds light on the non-thermal emis- 
sion processes (e.g. synchrotron and inverse-Compton pro- 
cesses), produced by a population of relativistic electrons 
in the jet. The degree of correlation between the spectral 
index a and the flux in various bands, through a least- 
square linear regression was checked, and values character- 
ized by large errors and bad x 2 were rejected. We found 
that the spectral index a of PKS 0735+178 varies be- 
tween 2.11 ± 0.26 to 0.66 ± 0.20, with an average value 
of 1.25 ± 0.15 in the last ten years. These values for a are 
ro ughly in agreement with the valu es calculated previously 
in lFiorucci. Ciprini. fc Tostil (|2004D using only the Perugia 
Observatory data set. 

In FigJS] (upper panel) the temporal behaviour of a is 
represented in comparison with the flux light curve in the 
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Fig. 5. Upper panel: temporal behaviour of the optical 
spectral index a and the i?-band flux in PKS 0735+178. 
This plot shows that the long-term variability is essen- 
tially achromatic, even if flares can imply spectral changes. 
Middle panel: dependence of the optical spectral index a 
on the flux intensity in the i?-band. The scatter in the data 
can be explained as produced by statistical fluctuations due 
to instrumental and calculation errors, and by the intrin- 
sic scatter of a during the more rapid and larger flares. A 
weak linear correlation is outlined. Lower panel: the scatter 
plot between the fluxes in R and B bands shows, as well, 
a linearly correlated emission (linear correlation coefficient 
r B -R = 0.90 ± 0.18 and slope 0.64 ± 0.04). 



better-sampled i?-band. The long-term variability seems 
essentially achromatic, and there is no obvious correlation 
between the light curve and the spectral index, whereas 
flares and short term variations can imply spe ctral changes. 
This i s the same behaviour found in BL Lac (IVillata et alJ 
20021 l2004h and in S5 0716+71 (|Ghisellini et all 11997k 
Raiteri et al.l feOOS). Unfortunately there is almost no spec- 



tral information during the 2001 outburst (lack of B and V 
data), to check a spectral flattening. The few data suggests 
a rather constant spectral index. In the same FigfS] (middle 
panel) the scatter plot between a and the flux is reported. 
Data dispersion is evident, and can be explained as statis- 
tical fluctuations due to uncertainties and to scattering in 
a during the more rapid and larger flares. Such linear cor- 
relation seems weak and a general spectral flattening is not 
detected clearly. Uncorrelated random fluctuations in the 
emitted flux might introduce a statisti cal bias, due to the 
spect ral index dependence by the flux dMassaro fc Trevesel 
Il996f ). but values computed for the central frequency (close 
to the i?-band, as plotted in Fig. [5J, can be considered unbi- 
ased and representative of the brightness state. The spectral 
index a showed consistent variations even when the light 
curve has a rather small variations. On the contrary in the 
lower panel of Fig. the scatter plot between the fluxes 
in the R and B bands shows a well correlated emission as 
expected (linear correlation coefficient rB-R = 0.90 + 0.18 
and slope 0.64 ± 0.04), without a detectable curvature. 

During well-defined and large flares at X-ray bands 
(especially observed in HBL), the X-ray spectral index 
versus the flux freque ntly displays a characteristic loop- 
like pattern (see, e.g . iGeorganopoulos fc Marscherl 119981 : 
iKataoka et al.1 



^ e.g . IL 
200$ [Ra 



avasio et al. 1 120041 ) . That "patterns 



outline a hysteresis cycle arising whenever the spectral slope 
is c ompletely controlled by radiative cooling processes (see, 
e.g. | Kirk. Rieger. fc Mastichiadisfl9 98: Bott cher fc Chiang 
2002). In fe w sources this feature was found in the optical 
regim e too (jFiorucci. Ciprini. fc Tostil 12004 ICiprini et al. I 
12004( 1. Consequently we can claim that around and beyond 
the synchrotron peak frequency, the behaviour of the LBL 
sources during flares in the optical band, is scaled in fre- 
quency but possibly very similar to the behaviour of the 
HBL in X-rays bands. 

In our 10-year light curve, PKS 0735+178 showed sev- 
eral moderate-amplitude outbursts, wider bump of longer 
duration, and a general flickering or shot-noise type of vari- 
ability on mid-term scales. The evolution of a as a function 
of the flux is erratic and did not show evident hysteresis 
loops caused by non-thermal cooling. In Fig.[S]thc evolution 
of a during 3 contiguous observing periods (from January 
29, to March 21, 2000) is reported as example. A rough 
loop-like behaviour is hinted, meaning that radiative cool- 
ing can dominate the optical SED also during mild-flaring 
activity. Consequently variations at higher frequency band 
could lead those at the lower frequency bands during both 
the increasing and decreasing brightness phases, reflecting 
differences in electron cooling times. 

The rather limited amplitude of the optical variability 
in the epochs of Fig. El the possible superimposition of dif- 
ferent emission processes in the optical band, the under- 
sampling and the error propagation in the a calculation, 
can be the main reasons for the lack of well defined loops 
in such a vs flux diagrams. Our data are not sufficient to 
make a final judgement, and an improved multi-band mon- 
itoring and a better data sampling would probably clarify 
the existence of that patterns also during mild variability 
in this object. 

8. Temporal variability analysis 

Time series analysis (evolved from both signal-processing 
engineering and mathematical statistics) provides very use- 
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Fig. 6. Evolution of the continuum optical spectrum of PKS 0735+178 (spectral index a) as a function of the flux in 
i?-band, during 3 contiguous periods of flickering variability with moderate flaring activity observed by the same telescope 
(Perugia Observatory, in period: January 29 - March 21, 2000). An approximative loop- like behaviour of a is hinted in 
the first (clockwise) and in the third (counterclockwise) patterns (see text for details). 



ful methods to study blazar variability. These methods al- 
lows to explore and extract temporal signatures, structures 
and characteristic timescales (the powerful scales of vari- 
ations), duty cycles (the fraction of time spent in an ac- 
tive state) and trends, to determine the dominant fluctua- 
tion modes and the power spectrum of the signal. Moreover 
time series analysis allows to detect and study auto/cross- 
correlations, time lags, transient events, periodicity and 
composite modulations, scaling and coherency, oscillations, 
beatings and instabilities, intermittence and drifts, dissi- 
pation, dumping, long-memory patterns and self-similarity, 
resonance and relaxation processes, random and determin- 
istic features, linear and non-linear processes, stationary 
and non-stationary activity, as well as to perform filter- 
ing and forecast. The analysis of the flux evolution over 
time in a blazar, joint with the multiwavelength and cross- 
correlation analysis, provides crucial information on the lo- 
cation, size, structure and dynamics of the emitting regions, 
and shed light on physical mechanisms of particle acceler- 
ation and radiation emission. 

In this section a quantitative analysis of the optical vari- 
ability observed in PKS 0735+178 is performed using 7 dif- 
ferent methods in 13 different light curves. The aim of this 
work is to examine in detail the optical behaviour on long- 
term timescales (months/years) using the whole 1906-2004 
historical light curve, the pre-1970 portion and the best 
sampled 1970-2004 part, and to explore mid-term scales 
(days/weeks in intervals < 200 days) using our 10-year 
monitoring data set with improved sampling. Each seasonal 
light curve in the best sampled i?-band is analyzed sepa- 
rately (a part of the first 2 pilot seasons, where we obtained 
only few observations). For the first time a rather com- 



prehensive temporal analysis of the optical variability was 
performed over 3 decades in time in this peculiar blazar, 
investigating scales between about 2 and 30 years regard- 
ing the historical dataset, and between few days and about 
100 days in the seasonal monitoring (the duration of season 
with data points spans between 144 and 203 days). 



T PKS 0735 + 178 t".4y 
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Fig. 7. Discrete Autocorrelation Function (DACF; data 
bin: 7 days, DACF bin: 150 days) and Lomb-Scargle 
Periodogram (LSP; dashed line indicates the threshold of 
false alarm probability [FAP] fixed at 1%) of the pre-1970 
light curve. A characteristic timescale around 11 years is 
hinted also by this poorly sampled portion of the historical 
light curve. 
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Fig. 8. Panels from left to right and following below: the whole historical (1906-2004) flux light curve of PKS 0735+178 
in S-band and plots from the related time-series analysis. First order Structure Function (SF; data bin: 7 day, SF bin: 
60 days) in log-log representation, Discrete Auto Correlation Function (DACF; data bin: 7 day, DACF bin: 90 days), 
Lomb-Scargle Periodogram (LSP; dashed line is the 1% false alarm probability [FAP] threshold), "Clean" implementation 
of the Discrete Fourier Transform (CDFT), Phase Dispersion Minimization function (PDM), plane contour plot of the 
wavelet scalogram (i.e. the two-dimensional energy density function \CWT(t,T)\ , CWT being the Continuous Wavelet 
Transform computed using a Morlet waveform, t and T the time and period scale respectively). In the last panel the 
periodogram of the synthetic light curve constructed upon the empty gaps using different window functions (GWFP; gap 
threshold 1 year). In the CWT scalogram the power spectral density is represented by filled-colour contour levels, while 
the thick black contours are the 90% confidence levels of true signal features against white/red noise background, and 
the cross-hatched regions represent the "cone of influence" , where edge effects become important. The descriptions of 
the plots above can be applied to the following FigEl FigfTUl FigfTll and FigfTJl Issues and results from these diagrams 
are described in the text. 



The following methods, optimized or adapted for un- 
evenly sampled timeseries, are applied: the first-order 
Structure Function (SF), the Discrete Auto Correlation 
Function (DACF), the Lomb-Scargle Periodogram (LSP), 
the Discrete Fourier Transform in the "Clean" implementa- 
tion (CDFT), the Phase Dispersion Minimization (PDM), 
the scalogram of the Continuous Wavelet Transform 
(CWT) using different waveforms, and the periodogram of 
the synthetic light curve constructed on the empty gaps, 
using different window functions (Gaps Window Function 
Periodogram GWFP). 



The SF is equivalent to the power spectral density func- 
tion (PSD) of the signal calculated in the time domain 
instead of frequency space, which makes it less depen- 
dent on sampling problems, like windowing and alias (see , 



e.g. Rutman 119781 : ISimonetti. Cordes. fc Heeschenl 119851 : 



ISmith et al.lll993T ). The first order SF represents a measure 
of the mean squared of the flux differences (Fj — F i+ At ) of 
N pairs with the same time separation At: 



1 N 

SFW(At) = -^(F i - J F i+At ) 



(1) 
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Fig. 9. Panels from left to right and following below: the whole historical (1970-2004) flux light curve of PKS 0735+178 
in -B-band and plots from the related time-series analysis. SF (data bin: 1 day, SF bin: 15 days), DACF (data bin: 1 
day, DACF bin: 60 days), LSP, CDFT, PDM, Morlet-CWT scalogram, and GWFP (gap threshold 0.3 years). Issues and 
results from these diagrams are described in the text. 



The general definition involves an ensemble average. Deep 
drops in the SF shape means a small variance and pro- 
vides the signature of a possible characteristic time scales, 
but a wiggling pattern and fake breaks (that can indicate 
false time scales) are common when the sampling is not 
sufficient. Typically, the SF increases with At in a log-log 
representation, showing an intermediate steep curve, whose 
slope b is related to the power law index a of the PSD by 
the relation a = 1+6 (a typical PSD has indeed a power-law 
dependence P(f) oc l/f a on the signal frequency / = l/t). 
The maximum correlation ti mescale is reached when the 
SF is constant for longer lags (|Hughes. Aller fc Allerlll992l : 
lLainela fc Valtaoidll993[ ). 



The DACF allows to study the level of auto-correlation 
in unevenly sampled data sets (se e, e.g. lEdelson fc Krolik! 
Il988t iHufnagel fc Bregmanl Il992t ) without any interpola- 
tion or addition of artificial data points. The pairs (Fi, Fj) 



of a discrete datasets are first combined in unbinned dis- 
crete correlations 

UDACFy = (^-<^»(^-<^» > (2) 

ffivfTj? 

where < F > is the average values of the sample and <jf, 
the standard deviation. Each of these correlations is asso- 
ciated with the pairwise lag Aty = tj — U and every value 
represents information about real points. The DACF is ob- 
tained by binning the UDACF^- for each time lag At, and 
averaging over the number M of pairs whose time lag Atij 
is inside At, i.e.: DACF(At) = l/M^-UDACFy. The 
choice of the bin size is governed by a trade-off between 
the desired accuracy in the mean calculation and the de- 
sired resolution in the description of the correlation curve. 
A preliminary time binning of data usually leads to bet- 
ter results. The number of the real points per time bin can 
vary greatly in the DACF, but data bins with an equal 
population can be built together with Montecarlo estima- 
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tions for peaks and uncertainties, as done in the Fishe r 
z-transformed DACF method (ZDACF. [Alexander Ifl997h . 

The LSP is a technique analogous to the Fourier analysis 
for discrete unevenly sampled data trains, useful to detect 
the strength of harmonic compo nents with a c ertain angular 



frequency to = 2ir f (see, e.g. Lombl 19761: Scargld 1982t 



iHorne fe Baliunaslll986t iPapadakis fc Lawrencdll993 ). 



In the CDFT method, first a "dirty" discrete Fourier 
transform (DFT) for unequally spaced data is calculated 
and then an inter active "cleaning" of the dirty DFT is per- 
forme d (see, e.g. lHogbornlll974lRoberts. Lehar. fe Dreherl 
Il987t iFosterl ll995|). The CDFT method is a complex 
and one-dimensional version of a deconvolution algorithm 
widely used in 2-dimensional image reconstruction. This 
technique provide a simple way to understand and remove 
false peak artifacts introduced by empty gaps. This method 
is effective especially in describing and recognizing multi- 
periodic signals. A standar d (uncleaned) D FT method was 
implemented previously bv iDeemingl (l!975h. 

The PDM method (see e.g lLafler fe Kinmanl 119651: 
lJurkevichllT97ll : IStellingwerill978h try to minimize the vari- 
ance of data at a constant phase with respect to the mean 
value of the light curve. If a trial period is close to a real 
period, the scattering of data against the derived mean in 
the light curve constructed on such phase (the light curve 
folded on such period) is small. The PDM method has no 
preference for a particular periodical shape, it incorporates 
all the data directly into the test statistic and it is well 
suited for small and randomly spaced samples. A value is 
statistically significant when the PDM drops towards zero. 

The Wavelet method is used to transform a signal into 
another representation able to s howing the inform a tion in 
a more useful shape (see, e.g . iDaubechiesI Il992t IFosterl 
119961: iPercival fc Waldenll2002t) . Wavelet transforms (WT) 
permits a local decomposition of the scaling behavior in 
time for each quantity (in contrast to the usual methods 
based on the Fourier analysis), allowing the signal features 
and the frequency of their "scales" to be determined si- 
multaneously. Hence it is a useful tool especially to detect 
typical timescales and identify signals with exotic spectral 
features, transient information content and non-stationary 
properties. WT are defined following the Fourier theory, 
but wavelets can be formally described as localized, oscil- 
latory functions whose properties are more attractive than 
sine and cosine functions. 

WT is computed at different times in the signal, us- 
ing mother wavelets (orthogonal base functions localized 
in both time and pulse spaces) of different frequency and 
convolved on each occasion. In this way the power spec- 
trum (i.e. the modulus of the transform value) on a two 
dimensional location-frequency plane is obtained (the so 
called wavelet "scalogram"). A continuous WT of a one- 
dimensional (ID) time-series is computed as a complex ar- 
ray at different times, the real component being the am- 
plitude and the imaginary component providing the phase. 
The square of the transformed modulus gives the wavelet 
power spectrum in function of both time and harmonic fre- 
quencies. In our analysis w e chose the Morlet complex- 
valued waveform (see, e.g. iFargd fl992T h composed of a 
plane wave modulated by a Gaussian envelope of unit 
width: 



Mr) = -n=e«*> T t 

V 



72 



(3) 



where r is the non-dimensional time parameter and u>o the 
non-dimensional frequency. Such continuous WT is convo- 
luted on the discrete sequence of the time series {Fi} with 
scaled and translated versions of tpo( T )- It is considerably 
faster to calculate such continuous WT in Fourier space: the 
convolution theorem allows to do all the N convolutions for 
a given scale simultaneously and efficiently in Fourier space 
N being the nu mber of points in the time series) usi ng a 
standard DFT (jKaiserlll994t IPercival fc Walden Il2002f) . 

The temporal analysis of the PKS 0735+178 optical 
dataset is performed using all the methods mentioned above 
separately on each light curve with results summarized in 
Table |3l In particular the diagrams from the time-series 
analysis of 5 light curves (the whole historical 1906-2004 
and best sampled 1970-2004 i?-band series, the seasonal 
IV, VII and X i?-band light curves from our dataset) are 
reported in Fig. |SJ FigHFigGl FigOU and FiglH In the 
analysis of mid-term scales, the problem of spurious arti- 
facts given by seasonal gaps with no data (solar conjunction 
with the source) was avoided studying each seasonal light 
curve separately. The GWFP of the best sampled 1970-2004 
light curve (Figl9] last panel, bottom, right) shows indeed 
only one powerful peak placed at exactly 1-year scale. This 
represents the periodical yearly recurrence of the seasonal 
gaps. 

The whole 1906-2004 light curve (Fig. G§ is patently af- 
fected by substantial differences in data sampling, by void 
gaps, by a poor sampling earlier than 1970 and a long 
empty interval (1958-1970). Nevertheless several signal fea- 
tures and characteristic timescales are pointed out by the 
different techniques (using both binned and unbinned se- 
ries): about 8.6y, 12-13y, 25y and 34y (see the summary 
reported in Tabl^3|). The 13. 7y timescale (pointed out by 
nice features in the LSP and CWT functions) is the same 
value cla imed as the major c ompo nent of a multi-p eriodical 
trend bv iQian fc faol (|2004f ) and lFan et all (|l997l ). on the 
other hand the 8.6y timescale (suggested by the DACF and 
LSP) is probably to be ascribed mainly by the temporal be- 
haviou r after 1970 (Fig(9j), and it is reported by lQian fc Tad 
(|2004D too. Longer duration scales (e.g. 34y) are difficult to 
be set out with confidence. Fake signal features, due to re- 
currences and temporal patterns given by the empty gaps, 
occurred at scales shorter than 8 years only (see the GWFP 
plot), therefore minor and fainter hallmarks in the SF, LSP, 
CWT corresponding to such shorter scales are neglected. 

The best sampled portion of the historical light curve 
(1970-2004, Fig|9j) spans 33 years and has a quite fair con- 
tinuous, regular and long-term coverage: 867 nights with 
1 datapoint at least, an average number of data points 
per night of 1.8, an average gaps among data of 7.8 days 
and a maximum gap of 1.6 years (Tablc[3|). The more rel- 
evant characteristic timescales suggested by the different 
temporal analysis methods are around 4.5y, 8.6y, and 12. 5y 
(other few values such 3.5y, 7.4y, 11. 8y years could be traced 
back to the previous mentioned, if we consider the finite- 
resolution accuracy of these methods). The 4.5, 8.6, 12.5 
years scales might be considered multiples harmonic sig- 
natures of a fundamental (coherent, absolute, or transient 
periodical, or again with drifting duration) component of 
about 4 years (see also the spline visual envelope in Fig[T2]). 
A characteristic time scales of about 4.8 year s was previ- 
ously claimed also bv I Webb et alJ (|l988l ) and lSmith et al.l 
(198 7]), while the 8.6 years scale was recently suggested 
by lQian fc Tad (|2004h . On the other hand no evidence for 
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Fig. 10. Panels from left to right and following below: the i?-band light curve of PKS 0735+178 in the IV observing 
season (Sept. 1995, Apr. 1996) obtained by our monitoring programme, and related functions produced by the time- 
series analysis. SF (data bin: 1 day, SF bin: 2 days), DACF (data bin: 1 day, DACF bin: 2 days), LSP, CDFT, PDM, 
Morlet-CWT scalogram, and GWFP (gap threshold 5 days). Issues and results from these diagrams are described in the 
text. General results of the analysis of all the seasons from our dataset are summarized in TableO 



th e periodica l signa ture of 14.2 years previously claimed 
bv lFan et alJ (119971 ) are observed, but only weak hints for 
scales in the 11.6-13.5 years range are found out. The power 
spectral density in the l// a regime, shows a slope index a 
between 1.5 and 2 (i.e. next to a pure shot-noise behaviour). 
The GWFP show a very powerful fake signature at 1.0 years 
as expected, produced by the recurrent 1-year gap between 
subsequent observing seasons. This artifact due to sampling 
it is not completely neglected by the methods used (see e.g. 
the residual peaks around uj ~ 6 i.e. ~ 27r in the LSP plot), 
therefore it is very useful to develop and make use of the 
GWFP technique in conjunction with the other time se- 
ries methods. No other (longer) fake features due to the 
irregular sampling are pointed out by the GWFP, there- 
fore all the other characteristic timescales claimed in this 
light curve can be considered due to real variability. In the 
CWT scalogram a significant and localized pulse in power 
is visible (with a scale around 4.8 years and located in cor- 
respondence of the big 1977 outburst). Another localized 
bump gives a scale of about 7.4y. An elongated and lower 



intensity band in the CWT scalogram, corresponding to 
the last epochs (about 1994-2004) indicates again a coher- 
ent time scale between the values 8.2-8.9 years. Results are 
quite similar using different suitable CWT mother func- 
tions as the Mexican-hat and Paul waveforms in this light 
curve. 

As example, in the time series analysis of the IV observ- 
ing season of our i?-band dataset (Sept. 1995, Apr. 1996, 
FigllOp a main modulating trend quite monotonically de- 
creasing is observed. This is stated for example by the shape 
of the SF in the logarithmic plot (second panel, top right 
of Fig fT0|) being quite linearly monotonic and without any 
turnover flattening (given by the reaching of a maximum 
correlation timescale). In this general trend 4 or 5 mod- 
erate and secondary oscillations can be visually identified 
(possibly related to the unique and weak signal of a char- 
acteristic timescale between 28-34 days, as hinted by the 
DACF, PDM, and CWT). Other characteristic scales are 
not displayed, and the artifact noise given by the irregular 
gaps, is important only at timescales below 13 days (see the 
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Fig. 11. Panels from left to right and following below: the i?-band light curve of PKS 0735+178 in the VII observing 
season (Oct. 1998, May 1999) obtained by our monitoring programme, and related functions produced by the time-series 
analysis. SF (data bin: 1 day, SF bin: 4 days), DACF (data bin: 1 day, DACF bin: 2 days), LSP, CDFT, PDM, Morlet- 
CWT scalogram, and GWFP (gap threshold 5 days). Issues and results from these diagrams are described in the text. 
General results of the analysis of all the seasons from our dataset are summarized in Tablcj3l 



GWFP last panel, bottom right). The power spectral den- 
sity in the l// a regime, shows a power index a = 1.97+0.25, 
i.e. a temporal variability mode like the shot noise (brown 
noise) signal. 

In the i?-band light curve during the VII observing sea- 
son (Oct. 1998, May 1999, Fig. ITT)) a brightening stage (be- 
tween about JD = 2451170 and JD = 2451265, i.e. 95 
days long) looking as produced by two blended main flares 
of about 50 days duration, might be supposed by a visual 
inference. Such values (about 95 and 50 days) are pointed 
out indeed by the SF, DACF, LSP, PDM and CWT func- 
tions (Tabled as characteristic timescales (scales where is 
more power in the signal) . The SF and PDM methods sug- 
gests also a possible timescale about 30 days. Fake features 
given by the irregular gaps are important only at timescales 
below 21 days (GWFP last panel, bottom right of Fig. ITT]) . 
The power spectral density function shows a power index 
a = 1.64+ 0.09, i.e. a temporal fluctuation mode placed at 
halfway between the flickering and the shot noise behaviour. 



The plots from the analysis of the X observing season 
data (Oct. 2001, May 2002, FigUJ) obtained by our mon- 
itoring programme, show an high state with two relevant 
flares, followed by a low and rather oscillating phase (af- 
ter date JD — 2452273). Unfortunately there was a long 
observing gap (about 25 days between JD = 2452218 and 
JD = 2452243) during the more active phase of this season, 
that is also the brightest optical state of the source recorded 
in our 1994-2004 database. Characteristic timescales of 
about 52-55 days are suggested by the LSP, CDFT, and 
PDM, while a 41 days scale is also found by the DACF and 
LSP (Tablc|3]) . Fake features given by the irregular gaps are 
influential only at timescales below 18 days (GWFP last 
panel, bottom right, Fig [T2"|) . The power spectral density 
function shows a power index a = 1.46 ± 0.17, i.e. again a 
temporal fluctuation mode placed at halfway between the 
flickering and the shot noise behaviour. 

The summary of the temporal statistics and analysis 
results is showed in Table[3] About statistics the following 
information is reported: the observing interval and its dura- 
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Fig. 12. Panels from left to right and following below: the i?-band light curve of PKS 0735+178 in the X observing 
season (Oct. 2001, May 2002) obtained by our monitoring programme, and related functions produced by the time-series 
analysis. SF (data bin: 1 day, SF bin: 4 days), DACF (data bin: 1 day, DACF bin: 2 days), LSP, CDFT, PDM, Morlet- 
CWT scalogram, and GWFP (gap threshold 5 days). Issues and results from these diagrams are described in the text. 
General results of the analysis of all the seasons from our dataset are summarized in Tablcj3l 



tion (98.1 years in total, and between 144 days and 203 days 
regarding our observing seasons); the number of the effec- 
tive observing nights with one data point at least (867 data 
points in the 1970-2004 light curve, and between 20 and 62 
in our dataset observing seasons); the average number of 
data points per observing night in the interval (spanning 
between 1 and 2.3, implying no data clustering); the aver- 
age separation between 2 successive data points (3 days on 
average on all our 10 observing seasons); the maximum sep- 
aration (maximum empty gap) between 2 successive data 
points (no more than 25 days in the worst sampled season) . 
In time intervals where the SF slope can be recognized in 
the log-log representation, we calculated its power index b 
trough a linear regression. About the time scries analysis 
the following quantities are reported (when possible): char- 
acteristic timescales calculated by deep drops in the SF, 
the power law index a = 1 + b of the PSD in the l// a 
regime calculated by the SF, characteristic timescales in- 
ferred from the SF turnover to the plateau produced by 
times longer than the maximum correlation lag, timescales 



estimated from power peaks in the DACF/ZDACF, in the 
LSP and in the CDFT, timescales indicated by deep drops 
in the PDM and by peaks in CWT scalogram. 

The criteria adopted in order to possibly avoid to quote 
fake features and artifacts due to the irregular sampling 
and gaps, or edge effects, in Tabled are the following: 1) 
only timescales shorter than 1/2 or 1/3 of the interval du- 
ration were considered; 2) only the more relevant signa- 
tures in each method are preliminary considered; 3) among 
these we discarded the features that are not indicated by 
more than one method on the same light curve portion, 
when the signature is not particularly strong; 4) if some 
timescales, among the remaining, is still matched by the 
synthetic GWFP, these are discarded too. Most of the ap- 
plied methods take well into account the power of artifacts 
given by the irregular dataset and recurrences of gaps, as 
showed by Fig J8l9l flOpi H and [12] (see the comparison of the 
first 7 panels with the last GWFP panel). However resid- 
ual spurious power can be still present in the functions (see 
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Fig. 13. Cubic spline interpolations of the best sampled 
1970-2004 historical light curve of PKS 0735+178. The con- 
tinuous (orange) spline curve is obtained with a data bin- 
ning of 600 days, while the other dotted/broken spline lines 
are obtained with binnings of 150 and 300 days. The inter- 
polations show 6 main peaks already visible with the 600- 
days bin. With this binning and considering the separation 
between the troughs of the pseudo-sinusoidal curve, we have 
cycles of about 6.0-6.5 years and 3.9-4.8 years. This raw vi- 
sual interpretation and spline envelope could be linked to 
one of the fundamental components (such as 4.4-4.8 years) 
found with the time series analysis (see Tab. [3]). Hence both 
the hypotheses of a multi-component recurrent trend, or a 
pseudo-periodicity (with period drifts) and modulating this 
light curve can be plausible. 



e.g. the spurious peaks of the LSP around u) ~ 6 i.e. T = 1 
year, in Fig[9]as mentioned). 

On long timescales the main temporal components we 
found can be grouped in three main ranges of values: be- 
tween 4.4 and 4.8 years, between 8.2 and 8.6 years and be- 
tween 10.8 and 13.2 years (Fig0Fig[51 Fig! and TabltiJ). 
These components could modulate the long term optical 
light curve of PKS 0735+178 with roughly cyclical os- 
cillations. A characteristic times cale of about 4.8 years 
was previousl y claim ed also by IWebb et al.1 (|l988f ) and 
ISmith & Nairl ll 1995ft whil e a 8.6 years scale was recently 

y\Qk 



suggested by iQian &: Tad (|2004ft . On the other hand we 
did not find any strong evidence of further and longer- 
duration (> 14 years) timescales re ported in literature (e.g. 
iFan et al.|[l997t IQian k. Taoll2004 and Sect©, but only 
weak hints of a 25 years and a 34 years signature. 

Some parts of the light curve could contribute with dif- 
ferent typical scales to the overall series, while data are 
treated in varying ways, with timescales suppressed or en- 
hanced as the relative weight of different segments of the 
light curve changes. The variability scale of 8.2-8.6 years 
could be an important finding and a real signature of a 
dominant and possibly quasi-periodical component anyway, 
because it was found in both the whole 1906-2004 series 
and its best-sampled portion. Moreover some uncertainty 
and statistical dispersion in the values found by different 
methods is expected, especially when the sampling is ir- 
regular, and here the dispersion range is small (0.4 years). 
In addition fast and/or isolated flares randomly occurred 
and uncorrelated to any general trend, can provide loud 
contributes to the power spectrum, disproving any period- 
icity hypothesis based on a long but under-sampled histor- 
ical light curve. The better sampled portion of the PKS 



0735+178 light curve did not disprove the characteristic 
timescale mentioned above, therefore it is reasonable to 
suppose a possible dominant period around 8.5 years. This 
hypotheses is open to future investigations based on pro- 
longed monitoring observations. In this view the shorter 
4.4-4.8 years scale found might be a submultiple of the pre- 
vious component. Hence this would be the real period of the 
fundamental harmonics (pointed out only by the best sam- 
pled 1970-2004 portion because of the sufficient sampling 
to detect it). A support corroboration of this conjecture is 
also provided by the spline interpolation reported in Fig |13l 
6 major maxima and cycles are outlined between 1970 and 
2004, with a duration between 3.9-4.8 years and 6.0-6.5 
years. This crude visual interpretation could be linked to 
the possible fundamental (and possible duration-drifting) 
modulating component of 4.5 years. Finally the group of 
longer-duration scales found in the broader range 10.8-13.2 
years are detected in each piece of the historical light curve 
(see Tab. [3]), but these values are probably too much scat- 
tered to mask a strict periodicity signature. 

On intervals shorter than 200 days, monitored by our 
observations in 10 subsequent seasons (from the III start- 
ing in Oct. 1994 to the XII ending in Feb. 2004) and ana- 
lyzed deeply with the methods mentioned above, there is 
no evidence for one single and pure periodical features, but 
there are signatures of several characteristic timescales of 
mid duration, commonly found in different observing sea- 
sons. These "recurrent" and "common" timescales are dis- 
tributed in few groups of values: 18 days, 24-25 days, 27-28 
days, between 40 and 42 days, between 50 and 56 days, 
65-66 days, between 76 and 79 days, and 95-96 days. In 
particular the timescales of 27-28 days are found in 3 ob- 
serving seasons (might be related to the synodical month 
interference), timescales between 50 and 56 days are ob- 
served in 6 seasonal light curves, and timescales between 
76 and 79 days are detected in 4 seasons. About these 
timescales, several hypotheses can be proposed. 1) These 
temporal signatures could be the result of a rough multi- 
periodical behaviour given by the superimposition of few 
harmonic components (spanning from about 2 dozen of 
days to about 100 days). 2) They could be produced by 
pseudo-periodical cycles, with a drift of the period duration 
around a fundamental value of 27 days for example (50-56 
and 76-79 days could be though as rough multiples in this 
case). 3) Such characteristic timescales could be produced 
by different periodical stages of transitory nature, with a 
sort of time-localized periodicity surviving only for limited 
epochs. 4) Again they could be the result of a variability 
mode endowed of few typical duty-cycles showing similar 
and recurrent peak shapes and durations, but occurring 
at random times. An improved and continuous monitoring, 
with an higher precision photometry and a better sampling, 
will provide more significant statistical results about the 
optical variability of PKS 0735+178 at these timescales. 



9. Summary and conclusions 

Blazars are one of the most exciting class of AGN, and the 
primary know extragalactic sources emitting high energy 
gamma-rays. Variability monitoring is an important effort 
in the study of blazars for several reasons (even if well sam- 
pled light curves are a very big challenge to be obtained at 
optical wavelengths normally). 



Table 3. Summary of temporal statistics and characteristic timescales found in the optical light curves of PKS 0735+178 (when possible) using 7 different methods. 
Data sets investigated: 3 historical £>-band flux light curves (the complete 1906-2004 curve, the best sampled 1970-2004 part, and the 1906-1958 portion), and 10 
separated i?-band flux light curves obtained in each observing season of our monitoring program (from the III season started in Oct. 1994 to the XII season ended 
in Feb. 2004). The following data are reported in the table columns (from left to right) for each light curve: (1) optical band and name of the light curve interval; 
(2) duration of the light curve; (3) number of the effective observing nights N on with one data point at least; (4) average number of data points per observing 
night < n >; (5) average separation between two successive data points; (6) maximum separation (maximum empty gap) between two successive data points; (7) 
power-law index a of the power spectral density (PSD) in the l// a regime (/ = 1/T), calculated in the time-domain through the first order structure function 
SF; (8) characteristic timescales Td r calculated by deep drops in the SF; (9) characteristic timescales T to inferred from the SF turnover to a long-lag plateau in 
the log-log representation; (10) timescales T pe estimated from peaks in the discrete auto-correlation function (DACF); (11) timescales T pe derived by peaks in the 
Lomb-Scargle periodogram LSP; (12) timescales T pe derived by peaks in the 'clean" implementation of the discrete Fourier transform (CDFT); (13) timescales 
Tdr estimated from deep drops in the phase dispersion minimization function PDM; (14) timescales pointed out by peaks in the 2-D contour plot of the wavelet 
scalogram (the two-dimensional energy density function provided by the continuous wavelet transform CWT computed using a Morlet mother waveform). The 
errors in the timescales are difficult to be estimated in general, and the significance of a scale is high when it is relevant according to a method, and when several 
methods point out similar values. Columns from (2) to (5) show that the sampling is fairly regular without high data clustering. The criteria adopted in order to 
avoid possibly the quotation of fake features and artifacts due to the irregular sampling in this table, and results displayed here, are described in the text of the 
paper. 
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[days] 


B 1906-2004+ 


98.1y 


989 


1.7 


20.7 


12.78y 




11.6y,25y 




8.6y,24.7y 


8.6y,13.2y,33.7y 


34y 


8.2,12.6y,15.2y 


13.7y 


B 1906-1958+ 


52y 


122 


1.4 


114 


8.95y 




12.3y, 18. 5y 




11. 4y 


5.7y,10.8y 




11.6y 


10.9y 


B 1970-2004+ 


33.3y 


867 


1.8 


7.8 


1.63y 


1.5,2.0 


4.4y,8.1y,11.8y 


0.6y,1.5y 


3.5y,4.5y,8.6y 


8.4y,12.5y 


4.5y 


8.2y,12.6y,15y 


4.8y,7.4y 


R III Oct.94-Apr95 


191 


43 


1.8 


2.5 


20.8 




79 




18 






18,78 


25 


R IV Sep95-Apr96 


203 


62 


1.4 


2.3 


12.9 


1.97 ±0.25 


39 




28 








34 


R V Oct96-Apr97 


178 


53 


1.6 


2.1 


17.9 


1.77 ±0.2 


50,79 


36 




50,77 




77 




R VI Oct97-Apr98 


189 


52 


1.3 


2.8 


15.9 




32,66 




68 






33,66,97 


95 


R VII Oct98-May99 


189 


51 


1.2 


3.0 


21.0 


1.64 ±0.09 


96 


31 


54,96 


53,102 




30,54 


48,96 


R VIII Nov99-Mar00 


144 


36 


1.4 


2.7 


22.2 


1.84 ±0.12 


83 


78 








25 




R IX OctOO-MarOl 


153 


20 


1.4 


5.6 


31.2 




27,56 




33 


28 


28 




40,76 


R X Oct01-May02 


201 


62 


2.3 


1.4 


24.8 


1.46 ±0.17 


69 


65 


41 


41,54 


52 


55,81 




R XI Nov02-Apr03 


144 


42 


1.1 


3.3 


24.9 


2.34 ±0.12 


24,60 




28 


28 


27 


55 


40 


R XII Sep03-Feb04 


148 


33 


1.0 


4.6 


21.0 




18,55 




20,42 






56 


52 



f Time scales followed by "y" are expressed in years. 
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1) Short-term observations and multiwavelength (MW) 
snapshots obtained during broad but limited-duration cam- 
paigns cannot resolve all the puzzling questions about 
blazars, while long-term monitoring allows to investigate 
the behaviour and evolution of the emitted flux on dif- 
ferent scales. 2) The knowledge about time variability is 
crucial like the spectral variability in constraint emission 
models, and the observed behaviour on long scales can be 
cross-correlated with the MW observations usually avail- 
able on short timescales. 3) Radio-optical monitoring could 
be considered a farsighted effort too: it enable to construct 
long-term records of variability for several sources, useful 
for future researches. 4) Even if most of blazars seems to ex- 
hibit an irregular, uncorrelated and unpredictable temporal 
behaviour, their optical light curve shapes appear to be not 
trivial (sometimes signatures of long-term memory, tempo- 
ral self similarity and intermittence are displayed) , whereas 
in few known cases periodical/quasi-periodical components 
cannot be ruled out. 5) In addition, flare triggers and tar- 
get of opportunity alerts for space observatories and large- 
size telescopes are usually based on a regular and con- 
stant monitoring. 6) Time series analysis of sparse data 
sets (like blazar light curves), is a challenging, interdisci- 
plinary subject, being developed and applied on a wide va- 
riety of present-day research topics outside astrophysics. 7) 
Our fairly novel investigation and results on mid-term opti- 
cal timescales (days, weeks), could be also compared to the 
analysis of blazar gamma-ray light curves that will be pro- 
vided, at the same scales, by the forthcoming Gamma-ray 
Large Area Space Telescope (GLAST). In fact this high- 
energy space observatory will be a large field-of-view and 
all-sky monitor for flares and variability, allowing to record 
flux variations on over timescales > 1 day on hundreds 
of 7-ray blazar-like sources. 8) Moreover worldwide inter- 
national collaborations and the participation of amateur 
and schools/universities optical telescopes are now possi- 
ble (thanks to the development of CCD photometry and 
automation technology) meaning a valuable link for educa- 
tion and public outreach. 

With this in mind, during a long-term and painstaking 
optical monitoring programme, we have obtained, collected 
and analyzed the largest amount of optical data in 4 colors 
ever published on the prominent blazar PKS 0735+178, 
thanks to the collaboration of 3 professional observato- 
ries (Perugia, Torino and Tuorla) and 1 amateur facility 
(Sabadell). Furthermore a new VRI photometric calibra- 
tion of 7 comparison stars in the field of this blazar is pre- 
sented (Tab. [T]), joint with the reconstruction and analysis 
of the whole historical light curve (spanning now from 1906 
to 2004). These optical data are rather unique with respect 
to continuity, sampling and duration for this source, and 
the associated data analysis enough comprehensive, despite 
of natural difficulties (weathcr/sccing conditions, seasonal 
gaps, technical problems or limited manpower). About 500 
nights of observations, collected in more than 10 years (pe- 
riod 1993-2004), and providing 1332 new BVRI final data 
points on PKS 0735+178 are reported and investigated, 
aiming to a quantitative statistical description of the data 
set, a characterization of the multi-band behaviour, and an 
investigation of variability over 3 decades in time, for the 
first time in this blazar. 

During the last 10 years, PKS 0735+178 continued to 
show rapid and large-amplitude optical variations typical 
of blazars, even if the source remained in a rather low or in- 



termediate brightness state (mag R > 14), showing a mild 
flaring activity. However starting from the end of 1997 the 
source showed a clear increase of the average brightness un- 
til 2001, when an active phase occurred. In this decennium 
typical variations of about 2 mag in less than half-year are 
observed joint with a general wiggling pattern produced by 
a superimposition or succession of flares, and modulated 
by a slower (possibly achromatic and oscillating) long-term 
trend. The quiescent and mild-activity reported from 1994 
to the second half of 2000 in the optical band, was also 
pointed out recently in radio bands, as a perio d of quiescent 
flux a ctivity and highly twisted jet geometry (|Agudo et alj 
2006). In the whole ~100-years history of PKS 0735+178 
five brightest outbursts and active phases were observed: 
the last occurred in the period Feb. 2001-Oct. 2001, and the 
brightest outburst was observed in May 1977 (when the 
source reached its historical optical maximum B ~ 13.9). 

The analysis of the continuum optical spectrum of this 
blazar suggests, as expected, a correlation between the 
fluxes in B and R bands, while the long-term variabil- 
ity of the spectral index a appear to be essentially achro- 
matic and independent by the wavelength (Fig. [5|). High- 
amplitude and isolated flares can imply correlated spectral 
changes (usually a flattening, i.e. bluer when brighter), but 
our data showed usually a rather erratic evolution of a as a 
function of the flux, and few or weak hints of non-thermal 
signatures (see e.g. Fig. At these mid-term timescales 
and without an increased sampling, it is reasonable to ex- 
pect that the superimposition of pure synchrotron optical 
fluctuations and emission peaks cannot be easily disentan- 
gled from slower variability patterns produced by different 
mechanisms. 

A summary of the quantitative temporal analysis per- 
formed in each single observing season of our i?-band light 
curves is reported in TabOH Intervals < 200 days are in- 
vestigated discovering quite common characteristic scales 
of variability falling especially into value ranges of 27-28 
days, 50-56 days and 76-79 days. These signatures are the 
stronger contribution to the power spectrum, possibly pro- 
duced by correlated flares with typical duty cycles emit- 
ted by some charge/discharge-like mechanisms. On the con- 
trary hand these typical timescales might represent the ef- 
fects of single or few powerful random events that are com- 
pletely uncorrelated and infrequent. On other words, even 
if a characteristic and intrinsic timescale is found, this does 
not mean necessarily a discovery of a dominant modulation 
or periodicity. In fact such evidence could well be produced 
by events of random or transitory nature, or be misrepre- 
sented by a combination of different underlying components 
affected by an insufficient observing sampling. 

Moreover the shot-noise (Brownian/brown noise) be- 
haviour, pointed-out by the values of the power-law index 
a of the PSD (computed in each season and falling be- 
tween 1.46 and 2.34, see TabJ3]), reflects the nature of the 
variations and can be linked to the findings cited above. 
Red and brown noise are termed usually as 1//° (power- 
law decline) fluctuations, meaning that the occurrence of a 
specific variation is inversely proportional to its strength. 
Brownian variability can be produced by a sequence of ran- 
dom pulses endowed of long-term memory, where indepen- 
dent/discrete events and parallel relaxation processes (like 
shocks and knots, electron density fluctuations, magnetic 
field turbulence or plasma instabilities) might generate a 
succession of mild optical flares and oscillations with sim- 
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ilar duty-cycles, on mid-term timescales. We remark that 
fluctuation in the statistical moments and parameters like 
the PSD and varian ce may be intrin sic in red/brown noise 
processes (see, e.g. IVio et al.l [2005^ 1 . even if our analysis 
is essentially phenomcnological and model independent. In 
addition when we analyzed separately different segments of 
the light curve and applied methods like the wavelets, we 
took already into account non-stationarity problems. 

About the historical behaviour of PKS 0735+178, we 
found 3 main characteristic timescales having an exten- 
sion of about 4.5 years, about 8.5 years (possibly signa- 
tures of the same 4.5 years fundamental component) and 
between about 11 and 13 years (TablelS]). These scales 
could be the result of an oscillating and achromatic trend, 
modulating (with a pseudo-periodic or multi-component 
course) the long-term variability. Pseudo-periodicity can 
imply drifts in duration and modulation, while a multi- 
component trend has several different scales contributing 
in the composite modulation. In particular in the first hy- 
pothesis the majority of the the long-term characteristic 
timescales found might be multiple signatures of a base 
component slightly drifting and varying around the value 
of 4.5 years. The visual inference based for example on the 
light curves of Fig. [Hand Fig. [T31 can supports this statis- 
tical finding. A rather "humped" or "multi-bumped" cycli- 
cal activity is evident, possibly meaning a bimodal course, 
defined by an alternation of active and quiescent stages. 
The limited temporal range of observations having a suffi- 
cient sampling (1970-2004) does not allow to understand if 
this pseudo-cyclical activity is a stable or transitory phe- 
nomena, and several hypotheses can explain the observed 
trend. A quasi/pseudo-periodical behaviour (with a funda- 
mental component drifting/oscillating around a value of 4.5 
years); a multi-component modulation (by possibly differ- 
ent correlated mechanisms); a mere random or transitory 
occurrence; a combination of the previous scenarios (for ex- 
ample a mixture of a multi-component trend with a quasi- 
periodical component). The achromatic behaviour reported 
in section is in agreement with a dynamical model im- 
plying slower variations of the base level flux and long-term 
modulations, resulting for by variations in the beaming fac- 
tor of the jet. 

This kind of optical course could be better corre- 
lated to the radio flux behaviour and the twisted (maybe 
precessing) jet of PKS 0735+178 observed in details by 
some years. This peculiar blazar (both radio and X-ray 
selected, and also a gamma-ray EGRET source) shows 
quite slow variations in the radio bands and a com- 
plex morphology displaying several moving components. 
In previous literature the optical and radio history of 
PKS 0735+178 have already suggested a possible peri- 
odical activity. Periodicity in blazars has been debated 
for more than 40 years and several models were devel- 
oped to explain this prospect (for example precessing or 
helical jets and supermassive binary black holes, see e.g. 
Lehto fc Valtonenlll996HSillanpaa et al.lll996t IVillata et alj 
1998t iRieeer fe Mannheim! I l200Ct IValtaoja et all l200oT 
Ostorero. Villata. fc Raiteril l2004h . However only in very 
few (and well publicized) cases there is still a sufficient ev- 
idence of cyclical outburst (like in OJ 287), and usually 
there is a general scepticism about widespread periodicity 
in quasars/blazars light curves till now. On the other hand 
the search for supermassive binary black holes in extra- 



galactic sources, should become a major research topic in 
the next years. 

As final consideration we point out that our 10- 
year observations probably mapped 2 distinct phases: a 
stage of low or intermediate optical luminosity (1994- 
2000) and a phase of mild flaring activity (2001). That 
dual and possibly cyclic scenario might well be confirmed 
by the behaviou r of the radio flux and structure during 
the same years jGabuzda et al.l 119941: iGomez et all 120011: 
lAgudo et al.l l2006). The optical flux is believed to be mainly 
originated in the very inner regions of the jet, even if 
the stronger variability could occur much farther from th e 
central engine than previously expected (jMarscherl [2005) . 
Hence it is reasonable to conceive a correlation between 
the optical and radio flux on long-term scales, and during 
the most active events. An improved continuous and longer 
optical monitoring of PKS 0735+178, and the comparison 
of the optical flare events with the ejection and evolution 
of the superluminal radio knots based on long-term data 
records, will allow to shed light on the physics of the jet 
flow and flare mechanisms in this interesting object. 
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